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H ' Abstract 



We consider a model of discretized 2d gravity interacting with Ising spins where 
phase boundaries are restricted to have minimal length and show analytically that 
the critical exponent 7 = 1/3 at the spin transition point. The model captures 
the numerically observed behavior of standard multiple Ising spins coupled to 2d 

gravity. 
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1 Introduction 

Despite the progress in our understanding of 2d gravity coupled to matter, both in 
the continuum description [l|, 0, |^ and as a statistical model of random triangu- 
lations @, ^, ^, we have still no understanding of the region which may be most 
interesting, i.e. matter fields with central charge c > 1 interacting with 2d gravity. 
One beautiful feature of 2d gravity coupled to conformal matter theories with c < 1 
is that the influence of matter on the geometry of the 2d universes only depends on 
the central charge. This is expressed by the KPZ formula |I|, ^ ^. Recent numerical 
simulations suggest rather surprisingly that the same is true at least for some range 
of c > 1 0. It is often conjectured that c > 1 is associated with 7 > 0, where 7 
denotes the string susceptibility. In this letter we describe a solvable model with 
non-trivial critical behavior for which 7 = 1/3 and which explicitly realizes a sce- 
nario for c > 1 advocated in |^ (see also |^). For other attempts to cross the c = 1 



barrier we refer to |10, 11 



2 The model 

A discretized model of 2d quantum gravity is defined by summing over the trian- 
gulations one can obtain by gluing together equilateral triangles along their links 
to form closed surfaces, and where one assigns the weight to each triangulation dic- 
tated by Regge calculus. In the following we will restrict ourselves to triangulations 
with spherical topology. This means that we can ignore the Einstein-Hilbert term 
in the action. The cosmological term is proportional to the number of triangles so 
the partition function is given by: 

where T denotes the class of triangulations of the sphere^, Ct a combinatoric factor 
present for triangulations of closed surfaces, Nt the number of triangles in the 
triangulation T and 7/v the subclass of T whose triangulations consist of N triangles. 
Since 

Y^ 1- = e^oAf ;v7o-3 ^^ ^ 0{l/N)) (2) 

the model has a critical point /iq at which the continuum limit has to be taken. It 
is known that 70 = —1/2. 



^ We find it convenient to consider only non-degenerate triangulations, i.e. each triangle in the 
surface has three different vertices. In terms of the dual 4>^ graph it means that wc exclude tadpole 
graphs. 



It is easy to couple matter fields to this discretized model. As an example one 
can couple Ising spins to the geometry by placing one spin o"j on each triangle i. 
The partition function is now given by 

z(/., /5) = E e-^'^ E 7^ E e^p if E(^^^. - 1) I (3) 

where J2{ij) denotes the summation over all pairs of neighboring triangles for a given 
triangulation T and J2{ai} the summation over the spin configurations on T. We 
can couple multiple Ising spins to the discretized model in the same way. For a fixed 
lattice such an extension would be trivial since the models would not interact, but 
here they interact via their backreaction on the geometry. 

The partion function (0) depends on two coupling constants. For each (3 there 
is a critical point fidP) (and a 7(/3)) where the continuum limit has to be taken. 
In this way the possible candidates^ for a continuum theory will be labelled by the 
spin coupling /3. In the case of a single Ising model one can explicitly solve the 
model |l^, and the result is as follows: There is a critical value of /?, Pc, above 
which there is spontaneous magnetization and below which the magnetization is 
zero. Away from (3c the backreaction of the Ising spins on the geometry is effectively 
zero (7(/3) = 7o), while 7(/3c) = —1/3, in agreement with the KPZ formula for 
a c = 1/2 conformal theory. On a regular lattice the spin transition is a second 
order transition and the corresponding conformal field theory has c = 1/2. On 
the dynamical lattice we see a corresponding transition at jSc, but the transition is 
modified by the coupling to gravity, and in addition the backreaction of matter has 
modified the gravity theory. 

For c > 1, i.e. for more than two Ising models coupled to gravity, no analyt- 
ical solution of the model is known. However, numerical simulations indicate the 
following: 

1) There is still a critical point f3c below which there is no magnetization and 
above which the system is magnetized. O, ITB|, IT7 



2) Above f3c the geometry seems to be that of pure gravity. Below (3^ there is 
a region where the situation is unclear, especially for a large number of Ising 
spins, but there seems to be a range of values of (3 for which the fractal structure 
of the surface is very pronounced. For sufficiently small (3 the geometry again 
is that of pure 2d gravity [|15], 0. 



■^The mere existence of a critical point in the disretized model does of course not ensure that 
there is an underlying continuum field theory. 



3) For a sufficiently large number of Ising spins 7(/3c) will be positive |T^, |TB|, |^. 

4) The geometry at (3c seems to be a function of central charge c of the matter 
system alone, at least for some range of c > 1. By this we mean that the 
effective distribution of geometries is the same no matter whether we consider 
2N Ising spins at their critical point or A^ gaussian fields or any other conformal 
field theory with central charge c = N, coupled to gravity . 

In view of point 2) above it seems reasonable to attempt a description of the 
model (for a large number of Ising copies) in a region around jSc and for large (3 in 
terms of an effective model where baby universes, i.e. parts of the surface connected 
to the rest by a small loop, are completely magnetized and where a transition should 
reflect the alignment of spins in different baby universes. Below we shall deflne such 
a model (see also |§]) which can be solved explicitly by very simple means and show 
that it, indeed, captures the features listed above except for the pure gravity region 
for (3 small. 

The model is obtained by restricting the summation in (|^) over Ising spin con- 
flgurations as follows: Given a conflguration {o"j} on a triangulation T the corre- 
sponding spin clusters consist of the maximal connected subsurfaces of T on whose 
triangles the spins are aligned and we shall require that the boundary components 
of all spin clusters are minimal, i.e. of length 2. In other words, we require all 
phase boundaries separating spins of opposite sign to be loops of length 2. We shall 
henceforth indicate summation over spin conflgurations with this property by Yl'tcn}- 
Spin clusters with a boundary of length 2 are clearly the excitations of the spins 
with lowest energy. Our summation J^'iaA represent an self-consistent iteration of 
such spin excitations. 

Thus the one-point (or one-loop) function G{fi, (3) is deflned as 

G(/.,/3) = E e-^^-E'exp f f E(^^^. - 1)) (4) 

TeTi {<xj V^ (ij) J 

where 7i denotes the class of triangulations whose boundary is a loop consisting 
of two (marked) links and for later convenience the spins on the two boundary 
triangles are flxed to one. Similarly, one deflnes ra-point functions that are essentially 
derivatives w.r.t. yU of the one-point function. 

We note that G(/i, /3) is well deflned and flnite in a region of the (/i, /?) -plane 
that contains the domain of deflnition for the full Ising model coupled to 2d gravity 
and is contained in the half plane /i > /iq, where /iq is the critical point of pure 2d 



gravity, whose one-point function is given by 

More specifically, it follows by standard arguments that there exists a critical curve 
(/x',(/5),/5)with 

fio < fi'm < f^cW) (6) 

such that G(/i, /3) is analytic in /i, (3 on the right of this curve. It is also easy to 
show that /ic(/5) — i> /io as /? — i> oo. 

The susceptibility x{l^i P) is defined as 

r)C 

x(/^,/?) = -^ (7) 

and the string susceptibility exponent 7(/3) is given by 

X(/^, P) = /(/^ - l^'ciP)^P) + / _ f(g))-y{i3) + ^^ss singular terms, (8) 

where / is an analytic function. 

The corresponding quantities in pure 2d gravity will be denoted by Xo(/^) and 
7o. In particular 70 is given by eq. (0) and, as mentioned above, equals —1/2. 

3 Critical behavior 

In order to determine the critical behavior of the model we establish a self-consistency 
equation as follows. Let T G T[ be a triangulation with a spin configuration con- 
tributing to G{fi,P). Each phase boundary (of length 2) separates a baby universe 
from the rest of the triangulation and clearly any two such baby universes are either 
disjoint or one contains the other. Thus we may in a unique way cut off maxi- 
mal baby universes bounded by phase boundaries and close up the corrresponding 
boundary loops of length 2 in the remaining part of T to obtain a triangulation 
T &Ti (with the same boundary as T) all of whose spins are aligned. Conversely, 
we obtain a T G Tl and a spin configuration {ai} contributing to (^(yU, /3) by starting 
with a T G 7^ with all spins aligned, cutting it open along a set of links and finally 
gluing on baby universes with appropriate spin configurations along the correspond- 
ing loops of length 2 such that the spin on the two boundary triangles are oppositely 
oriented to those of T. 

Using this decomposition procedure and summing first over baby universes with 
spin configurations we obtain 

G(/x, /?) = E ^''""^ (1 + e-'^G(/^, /3))^^"' (9) 

TeTi 



where Lf denotes the number of hnks in T, i.e. Lf = 1 + \Nf. In eq. (^) the factor 
g-2/3 represents the couphng of a baby universe across the phase boundary to the 
rest of the surface and the factor 1 in the parentheses originates from the empty 
baby universe. We can now reexpress eq. (^ as0 

G(/x, /3) = E e-^"^" = G'o(/x), /x = /i - I log (l + e~'^G{^^, /?)) . (10) 

Note that the last equation can be written as 

/i = ^+^log(l + e-2'3Go(/x)) (11) 

which expresses /i in terms of known functions of /i and (3 since pure gravity can be 
solved. 

From eq. (p!OD and ([TT|) we get 

x(/x,/3) = Xo(/i)|^ (12) 

dl ^ e'^ + Gom) 

dfi e2/^-(|xo(/i)-Go(/x))- ^ ^ 

Since the string susceptibility exponent 70 = —1/2 < in the case of pure gravity 
both Go{fio) and Xo(Aio) are finite. This implies that there exists a Pc, given by 

e'^= = Ixoif^o) - Goifio) (14) 

such that the denominator in (|1^) is different from zero for all /i > fi'^iP) provided 
13 > (3c ■ By differentiating x(/i, (3) n times with respect to ^ we know from (P) that 
X*-" (/U, (3) will be singular for /i — ;> /ic(/5) provided n is sufficiently large. Performing 
the same differentiation on the rhs of (P^), using again (|T^, we see that the only 
chance for a singular behavior is that /i(/i'^(/5) = /io and in this case the leading 
singularity has to come from Xo {p)- We conclude: 

/i(/i^(/3)) = /io and 7(/5)=7o for /3 > /^o- (15) 

This is the phase where the model is magnetized, where the spin fluctuations are 
small and where the geometry of the surfaces is not affected by the spins. 



The number /?c (given by ([1^)) is characterized by being the largest j3 for which 
dn/dfl equals zero for some /2 > /xq, i-e. (by ([T3|) ) the largest /5 for which the 
equation 

e'' = |xo(/i) - G'o(/i) (= E (3A^t/2 - l)e-^^-) (16) 



4 



Strictly speaking we have added a factor (1 + e^^'^G(/i,/3)) on the rhs of (|^) since Lt ~ 2 — 
|jVy — 1. However, aU arguments presented in the foUows are vahd even if we did not include this 
factor. We have chosen to add it in order to simplify the formulas. 



has a solution for fi > fiQ. If we define fidP) by 

= (17) 



dfj, 

dfi 



Ac(/3) 

for P < Pc, then fldP) obviously solves (ITH) and we have 



fj^ciPc) = /^o, Aic(/5) > yWo for /3 < pc- 



Let us now assume that P < Pc- If we use eq. (|^) and (|T3|) this implies that 
xifJ^, P) will be singular for /2 ^ /ic(/?) due to the vanishing denominator on the Ihs 
of eq. (|1^). In fact, since fidP) > /Wo both Xoifi) and Go{fi) will be regular around 
ficiP) and we can Taylor expand the Ihs of (p!2D: 

/^ - /^c(/5) y^ - ^;^(/?) 



To derive the last equation we have used ([T^) and (|T^ which tell us that 

/x = fi'^iP) + const{fi - fic{P)Y + ■■■■ (20) 

We conclude that 7(/?) = 1/2 for /? < Pc- In this phase baby universes are dominant. 
Effectively we have branched polymers and the total magnetization of the system is 



zero 



Let us finally consider the system at the critical point Pc- This point is char- 
acterized by the fact that ftdPc) coincides with fiQ. Although the singularity of 
xifJ^, Pc) for /^ —* f^'ciPc) is still dominated by the zero of dfi/dfi we can no longer 
Taylor expand fi{fi,, Pc) around fidPc) (= /Wo) since the functions in (|TE|) are singular 
in fiQ. On the other hand we can use the known singular behavior of the pure gravity 
functions Gq and Xo at fio to deduce from (|13|), remembering that 70 < 0, 



^ ~ (/i - /i.(/3.))-^°, i.e. /x - fi'ciPc) ~ (/i - M(3c))'''^'- (21) 

From (jl^) we finally get, using (|2T|) 

c c 

^^^' ^"^ ^ jfi^fiMr^ ^ (/i-/^^(/?c))-^o/(-^o+^)' ^^^^ 

and we have derived the remarkable relation: 

7(/5c) = -^^ = \. (23) 

-70 + 1 3 



4 Conclusion 

We have seen above that our toy model has two generic phases, a magnetized phase 
for large /3 where the distribution of geometries coincides with that of pure gravity, as 
one naively would expect for large jS, and a phase where 7(/3) = 1/2. The exponent 
7 = 1/2 is that of branched polymers and the total magnetization of such a system 
is zero due to the linear structure of branched polymers [jl5|. At the critical point 
f3c for magnetization we have 7(/?c) = 1/3. 

Our model shows how the interaction between matter and geometry can lead to 
a string susceptibility 7 > 0. It agrees remarkably well qualitatively with numerical 
simulation of multiple Ising models except for a region of small f3 which shrinks 
with increasing multiplicity of the Ising models. The model is closely related to 
the matrix models studied in [ll^ and to the c —>■ 00 limit of multiple Ising models 



studied in |2y] . However, our approach has the virtue of being simple and avoids any 
use of matrix models by working directly with the spin excitations on the surfaces. 
Moreover, it highlights the general nature of eqs. (|12D and (0). 

It is a natural assumption that the full multiple spin model will have logarithmic 
corrections to 7 = 1/3 (if that is indeed the correct exponent) and that these will 
increase as the central charge c decreases from infinity. There is even some numerical 
support for this conjecture 0. 

T.J. acknowledge the support from the Fulbright Foundation and would like to 
express his thanks to for the hospitality at ITP and the Department of Mathematics 
at UCSB. 
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